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The coexistence of infinitely many attractors is called extreme multistability in dynamical systems. 

In coupled systems, this phenomenon is closely related to partial synchrony and characterized by 
the emergence of a conserved quantity. We propose a general design of coupling that leads to partial 
synchronization, which may be a partial complete synchronization or partial antisynchronization 
and even a mixed state of complete synchronization and antisynchronization in two coupled systems 
and, thereby reveal the emergence of extreme multistability. The proposed design of coupling has 
wider options and allows amplification or attenuation of the amplitude of the attractors whenever it 
is necessary. We demonstrate that this phenomenon is robust to parameter mismatch of the coupled 
oscillators. 

PACS numbers: 05.45.-a, 05.45.Xt, 05.45.Gg 

Multistability-the coexistence of many attractors is an intrinsic property of many dy¬ 
namical systems in physics, chemistry, biology, and in ecosystems and synthetic genetic 
networks. It has also been demonstrated in laboratory experiments, such as lasers and elec¬ 
tronics. Multistability reveals a rich diversity of stable states of a dynamical system and, as a 
consequence, such systems offer a great flexibility since they allow switching from one stable 
state to another either by the influence of noise or other intrinsic control mechanism. This 
is, particularly, important for the survival of species in ecosystem. Hence, understanding 
multistability and its control is an important issue. On the other hand, generating multista¬ 
bility is a challenge. We consider systems possessing not only a finite number of coexisting 
attractors but infinitely many of them. Such systems can be constructed by connecting 
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two dynamical systems with an appropriate design of coupling even if the isolated systems 
do not show multistability. We showed earlier that a partial complete synchronization is 
an essential condition for extreme multistability when a conserved quantity emerges. Be¬ 
sides complete synchronization, in physical, chemical and biological systems, other forms of 
synchrony such as antisynchonization or mixed synchronization may exist. Here we estab¬ 
lish that extreme multistability can be realized with partial antisynchronization, and even 
under the condition of mixed synchronization in two coupled systems, if we appropriately 
generalize the strategy of the coupling design. Additionally, we demonstrate that we can, 
in fact, increase and decrease the amplitude of all attractors in a desired way, whenever it 
is necessary. Furthermore, we establish the robustness of this phenomenon with respect to 
parameter mismatch of the coupled systems. 


I. INTRODUCTION 


f- .1] 


The coexistence of a multitude of different stable states is observed in many systems (cf. |l| for a review): bistability 
in neurocortical oscillations 2, 3], several possible overturning states in the thermohaline ocean circulation (THC) 
, multiple equilibria in a global ocean circulation |5J, different community compositions of species in ecosystems [61 
and various gene expression in synthetic genetic networks |7|. Multistablity is also seen in experiments in lasers and 


opto-electronic devices _8j, |9|, in condensed matter physics 


10l | and electronic circuits If 1 . From the theoretical point 


of view one can distinguish several system classes which possess coexisting attractors: weakly dissipative systems 
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12j, coupled systems 


13] or time-delayed feedback systems 


14]. Particularly in coupled systems the coexistence of 


multiple attractors is closely related to the coexistence of synchronized and desynchronized states. Multistability in 
coupled oscillators can be obtained (i) in cases where the single isolated oscillator already shows multistability as in 
two coupled rotors 


m- 


15] as well as (ii) in the case when the single oscillator possesses only one attractor as in coupled 


phase oscillators |16|. On the other hand, creating a large number of coexisting attractors is an interesting issue 
to investigate and its control is an even more challenging task. It has been shown that one can get an arbitrarily 
large number of coexisting attractors in weakly dissipative systems by choosing a very small dissipation 12]. Another 
possibility is to consider a drive response system where the driver is a conservative system, i.e. for each initial 
condition of the driver, the driving force on the response system is different leading to even an infinity of attractors 


If Alternatively, infinitely many attractors can be achieved when nonlinear systems are coupled in a specific way 
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to show partial synchrony 


isml- 


Theoretically and analytically^ 


multistability has been investigated in chemical reactions like the oscillatory 


22]. Even more complex multistable behaviours appears in chlorite-thiosulfate re- 


Belousov-Zhabotinskii reaction 
action 23]. Despite special care taken to ensure the same experimental conditions, this reaction exhibits different 


outcomes. One possible exp lanation of the experimental observation of seemingly infinitely many attractors has been 


offered by Wang et al. 


24( using the autocatalator model, where the inclusion of a buffer state creates the coexis¬ 


tence of different states. Later, Sun et al 


18] reported the existence of an infinite number of coexisting attractors, 


called extreme multistability (EM) in two coupled Lorenz systems by using a particular type of coupling. How to 
design a coupling scheme leading to extreme multistability in coupled oscillators, in a straightforward way, has been 
demonstrated in H and later in 


21]. The emergence of extreme multistability has also been evidenced recently in 


25| . The main characteristics of extreme multistability 


an electronic experiment with two coupled Rossler oscillators 

has been studied in detail in Q . This study reveals that two properties are essential for the appearance of extreme 


multistability in two coupled n-dimensional nonlinear systems: (i) the complete synchronization (CS) of n — 1 com¬ 
ponents of the two systems and (ii) the emergence of a conserved quantity C in the long-term limit t —> oo. Note 
that in CS of two coupled systems, where all the n components synchronize, the long term dynamics lies on a unique 


synchronization manifold 


26-28]. By contrast, in case of EM, only n — 1 components synchronize completely, while 


the differences between the n — th components of the two systems remain in a constant distance to each other. This 
constant distance C defines a synchronization manifold in which the dynamics takes place. This conserved quantity 
C, which appears only in the long-term limit t —> oo, depends on the initial conditions in a possibly complicated way. 
Since the value of C can take any real number, the whole state space is foliated into infinitely many synchronization 
manifolds with at least one attractor existing in each manifold. 

It is important to note, that the attractors in the synchronization manifold do not correspond to the usual definition 
of an attractor even in the Milnor sense Q . According to Milnor, an attractor possesses a basin of attraction 
of positive Lebegue measure. As pointed out above, the whole state space is foliated into infinitely many such 
hypersurfaces as time goes to infinity. This can be considered as a foliation of the state space into infinitely many 
leafs. Within the each hypersurface or leaf we obtain a usual attractor as in any other dynamical system but only 
with respect to this leaf. In the direction perpendicular to the leaf we have marginal stability. That means that 
the invariant set in the leaf is in fact a relative attractor [321 ] since it is stable in all directions within the leaf or 


hypersurface but marginally stable in the direction perpendicular to it. Since the basin of attraction of each attractor 
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consists in the simplest case of a 2n — 1 dimensional hypersurface, its Lebegue measure is zero in the 2 n dimensional 
state space, but the attractor would have positive measure with respect to the hypersurface or leaf when considered 
as a relative attractor. To avoid confusion, we would like to stick with the notion of an attractor instead of a relative 
attractor since this phenomenon of extreme multistability has been discussed earlier in several papers (cf. 


18 


m 


and we would like to refer to those papers using their notation. Throughout the whole paper all attractors should be 
understood as relative attractors. 

The general rule for defining the coupling function that successfully creates EM in two oscillators 


201 ] is based on 


Lyapunov function stability (LFS) 3.3]. Since different types of synchronization like antisynchronization or mixed 
type of sychronization may appear in coupled nonlinear systems, we address here the question, if EM can exist for 
partial antisynchronization (AS) or in a state of mixed synchronization (MS) besides partial CS. In a state of AS, 
all the pairs of state variables of the two oscillators are correlated by amplitude but are exactly opposite in phase 
(7r) while in the MS state, some of the pairs of variables are in CS and others in AS state. We make an appropriate 
generalization of the definition of the coupling function to achieve this goal and show that the number of choices of 
the coupling is not restricted but widely open. 

An interesting benefit of our specific coupling for two oscillators is a possible amplification or attenuation 33] of the 


amplitudes of the state variables of the coupled system. This option may help amplifying the smaller amplitude of 
the attractors noticeable or, on the contrary, limit an almost unbounded growth of an attractor by attenuating it. We 
make use of this possibility to change the amplitude of the emerging attractors in the state of EM. To demonstrate 
these novel features of EM we use two coupled Rossler oscillators as a paradigmatic model. 

So far, EM has only been demonstrated in coupled identical oscillators which are difficult to realize in experiment 
as shown in [25]. To provide a theoretical understanding of an experimental setup one has to consider slightly 


different oscillators and to show that EM persists under parameter mismatch. To this end, we test the robustness of 
the phenomenon of extreme multistability with respect to mismatch of parameters. 

The paper is organized as follows: In Sec. [IT] we describe the coupling scheme to achieve EM and extend it to partial 
AS and MS states as well as scaling the amplitude of attractors. We illustrate our method of creating EM with 
partial AS using numerical simulation of the Rossler oscillator in Sec. Mil Section IIVI deals with the effect of the 
parameter mismatch and finally we discuss our results in Sec. [V]. 











5 


II. COUPLING DESIGN AND CHOICE OF CONTROLLERS 


f. l 


We recall briefly the main ideas of the systematic design of coupling to achieve EM (for details cf. J20|) and extend 
it to partial AS and MS. The basic condition for the emergence of EM is: at least one pair of variables maintains a 
constant distance that emerges as a conserved quantity which is sensitive to changes in the initial state. To obtain 
partial AS we require all other variables to be antisynchronized, while for MS some pairs of variables should be in 
AS and the remaining ones should be in CS. For this purpose of generalization, we now introduce a scalar matrix a 
d in the definition of the error vectors of the coupled system. We start with two identical oscillators as described 
by x = F(x);x = [xi\i = 1,2..., n] and y = F(y),y =[yf,i = l,2...,n], where F : R n — >• R n . Assume the synchronized 
dynamics taking place on a manifold defined by x = ccy and the deviation from the synchronization manifold is 
described by the error e = x — ay. The constant matrix a (n x n) implements scaling up or down the comparative 


amplitude of the pairs of the variables of the coupled system 


33] and thereby takes care of increasing or decreasing the 


amplitude of the attractors. Additionally, it takes care of the type of synchronization. If the matrix a is the identity 
matrix then all the variables will be completely synchronized. However, all the elements (i,j = 1 n) of the 
a-matrix may not be positive natural numbers. They may take any real value (a^ £ R 1 ). For simplicity we focus 
here only on the matrices a which are diagonal. If the diagonal elements of the a-matrix are all negative but the 
values are -1, then the variables are in an AS state. The MS state can be created if some of the diagonal elements of 
the a-matrix are chosen to be +1 and some others are -1. Choosing the elements larger or smaller than unity, allows 
amplification or attenuation of the amplitude of the attractors of one system compared to another when the systems 
are coupled. The final form of the error dynamics is described as: e = G(x,y, a) = F(x) — aF(y), G : R n — > R n . 
To fulfill the above mentioned conditions, the error dynamics has to obey a specific targeted form which we denote by 
G(x,y,a). This desired error dynamics G(x,y,a) is realized by designing a set of controllers Ui(x,y) and u 2 (x,y) 
for coupling the two dynamical systems in such a way that it originates EM in the coupled system. Now the desired 
error dynamics e = x — ay = G(x, y, a) implies that G(x, y, a) — F(x) + aF(y) = Ui(x, y) — au 2 (x, y). According 
to a particular choice of the matrix a, one can choose the controllers Ui and u 2 in such a way that CS is obtained in 
p-pairs of variables and AS in the (n-l-p)- pairs of variables and, the remaining one pair maintains a constant distance. 
In this way any possible state, either a partial MS or a partial CS or AS can be realized. This generalizes the concept 


introduced in 


20l | to different kinds of synchronization and additionally, introduces an option to amplify or attenuate 


the amplitude of the attractors. 
As already mentioned in n 


the choice of controllers is not unique, but very flexible. Here we quantify the 
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flexibility of the choice of controllers. To derive the minimum number of different controller choices when coupling 
two n-dimensional oscillators, we start with 3-dimensional systems. First of all we mention that it is necessary to 
have a mutual coupling to achieve EM. For this reason we cannot set all the controllers of one oscillator to zero, i.e. 
uu = 0 or U 2 i = 0 Vi is forbidden. This leaves us with the possible choices (1) u \2 = 0; U 21 = 0 and M 23 = 0, (2) 
Mn = 0; U 22 = 0 and U 23 = 0 and (3) M 13 = 0; M 21 = 0 and U 22 = 0, where one controller of the first oscillator 
and two controllers of the second oscillator are set to zero. Alternatively we can also nullify two controllers of the 
first oscillator and one controller of the second oscillator which yields another 3 choices. Thus we can construct 
3 + 3 = 6 different types of controllers for one particular choice of error dynamics in 3D systems. Next we consider 
a 4-dimensional system for which we obtain ((f) + ( 2 ) + ( 3 ) = 14 choices of controllers from a simple counting of the 
possible choices for one specific error dynamics. In a similar manner, we determine the total number of choices of 
controllers for couplig two n-dimensional systems are (™) + Q) + (3) + + („"]} = 2" — 2. This follows from the 

relation (1 + x) n = 1 + (");r + + (ff)# 3 + + x n with setting x = 1. We emphasize that this minimal number 

of controllers is available for one particular choice of the error dynamics. 

However we can also choose U\ = p[G(x, y, a) - F(x) + aF(y )] and «2 = — (1 — p)[G(x, y, a) — F(x) + aF(y)] with 
p being any rational number between 0 and 1. As a consequence, in principle, infinitely many choices of controllers 
are available for coupling two systems to realize EM. This large coupling options opens up possibilities of a physical 
realization of the EM. 


III. EXTREME MULTISTABILITY: NUMERICAL EXAMPLES 

To illustrate the emergence of EM in a partial AS state, we consider two identical Rossler systems coupled through 
bi-directional controllers given by 


X1 

= ~X 2 - X3 + Mil, 

(la) 

X2 

= X\+aX2+U\2 1 

(lb) 

&3 

= b + x 3 (x 1 - c) +M13, 

(lc) 

m 

= -J /2 - 2/3 ++ 21 , 

(Id) 

m 

= 2/1 + °2/2 + U22, 

(le) 

h 

= b + y 3 (y 1 -c) +m 23 - 

(If) 
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(a) (b) 

FIG. 1: Multiattractor dynamics in coupled Rossler systems: 2D projection of different trajectories converging to different 
attractors for different initial states, xoi- (a) x\ vs xi and (b) yi vs 1 / 2 - an = —2.0 ,022 = —2.0 and 033 = —2.0. For 
comparison, two figures are drawn in the same scale. It confirms that one set of attractors is an amplified version of the other 
set and they are in antiphase. Other initial conditions are fixed at X02 = 0.0; X03 = 0.2; 3/01 = — 0.1; 2/02 = — 0 . 1 ;j/o 3 = 0.0; *01 
has only been changed from -4 to 4. 

mj (i=l, 2 and j=l, 2, 3) are controllers which are to be designed to define the coupling. The parameters of the 
uncoupled Rossler systems are set to such values that a chaotic attractor (a = 0.2; b = 0.2; c = 5.7) exists; no 
multistability occurs. We can now derive the coupling function by an appropriate choice of the controllers to realize 
EM with a partial AS or MS along with mutual amplification or attenuation of the coexisting attractors. We have 
already introduced the scaling matrix in the definition of the errors in the previous section. The errors have been 
defined by ei = Xi — auyi where a is a diagonal matrix (for simplicity) and its elements au(i = 1,2,3) are real 
numbers. Once a desired synchronization regime is targeted by a proper selection of an, we concentrate on the choice 
of controllers for realizing EM. We design the error dynamics in such a way that the coupled system is partially 
synchronized in any of the desired synchronization regimes (CS, AS or MS) in which, at least, one pair of variables 










evolves into a conserved quantity. For the present example we choose the desired error dynamics as 


ei = ~ 63 + 7/162, 

= ^263, (2) 

e3 = —ce3 — xie-2. . 

For a choice of Lyapunov function V = e| + e§, ei and e 3 tend to zero in the long run while ei emerges as a constant 
C depending on the initial conditions. Note that this choice of the error dynamics and the Lyapunov function is not 
unique as disucssed above. We can always make other choices provided they satisfy the basic conditions of partial 
synchrony and introduces an emergent conserved quantity depending upon the initial state. Using the difference 
equations ofWLaMVMbMeimW), and the error dynamics (J2]) we derive the controllers as follows 

Mil = 1/1 e 2 + X 2 — ot\y 2 + x 3 — aiy 3 — e 3 , 

U 12 = 0, 

M 13 = -x 2 e 2 - b + a 3 b - xix 3 + a 3 y 3 y 3 , (3) 

«2i = 0, 

x 2 e 3 — xi + a 2 y\ — ae 2 

u 22 — -, 

— 012 

M23 = 0 . 

Using these choices of the controllers, we obtain the coupled Rossler oscillators that exhibit EM, 

xi = -x 2 - x 3 + yi(x 2 - a 2 y 2 ) + x 2 - oi y 2 
+£3 - otiy 3 -x 3 - a 3 y 3 , 
x 2 = x\ + ax 2 , 

x 3 = b + x 3 (xi - c) - x 2 (x 2 - a 2 y 2 ) — b + a 3 b - x 3 x 3 

+ 031 / 12 / 3 , (4) 

Vi = 2/2 - 2/3, 

V 2 = 2/i + ay 2 

, * 2(^3 - o 3 2 / 3 ) - xi + o 2 2 /i - a(x 2 - o 2 2/ 2 ) 

I 

2/3 = b + x 3 (x\ - c). 


— 02 
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(a) 




(b) 


FIG. 2: (a) Maxima values of 22 as a function of initial condition £ 01 . The rest initial states are fixed as said before, (b) An 

emerging constant as a function of same initial condition 201 - 

The errors e 2 and e 3 go to zero in the long run since V = — ce\, e\ maintains a constant distance. For a choice 
of an = — 2 . 0 , a^i = — 2.0 and £*33 = — 2 . 0 , all the states are in antiphase and one set of attractors is expected 
to be larger by twice the other set of attractors. The conserved quantity still appears as a sum of two variables 
ei = X\ + 2yi = e\ (ef is the conserved quantity). 

Figures [Ija)-(b) show the attractors (aq vs 22 and yi vs yi plots) of the coupled system with variation of one initial 
state 2 oi- They clearly reveal the coexistence of a large number of attractors which are in AS state. The mutual 
amplification of the attractors (one is twice as larger as the other) is also clearly seen by comparison of the Figs. [U a ) 
and 02 b). 

The infinite number of coexisting attractors is in fact confirmed by the bifurcation-like diagram in Fig. by 
plotting the maxima or peak values of 22 for changing initial condition 201 . Many dynamical states (limit cycle or 
chaotic) appear for each initial state 201 . The emerging constant (xi(t en d)+2yi{t en d)) is also plotted versus 201 in 
Fig. [H(b). Here t en d is the final integration time which is chosen to be t en d = 500000. This confirms that the emergent 
constant follows a nonlinear relation with the initial state and thus it is challenging to predict the emergent attractor 
for a chosen set of initial conditions. We do not elaborate the MS state here with an example, since it is clear that, 
by using a different set of an, one may ensure MS which does not change the main algorithm for realization of the 


EM. 
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IV. ROBUSTNESS OF EM TO PARAMETER MISMATCH 


Finally we check the effect of parameter mismatch on EM. Since no two oscillators, in the real world, are exactly 
identical, there exists a mismatch in their parameters. This poses a problem for the experimental realization of EM 
by coupling two oscillators with a mismatch which may finally lead to the disappearance of EM. Recently, EM has 
been evidenced in a laboratory experiment [25] using two almost identical electronic analogs of Rossler oscillators 


employing a simple error dynamics 


20|. It was found that, as expected, the dynamics drifts continuously from one 


attractor to another qualitatively different attractor. This ’’drifting” of attractors appears due to an instability created 
by the parameter mismatch in the coupled oscillators. It was observed that if the coupled oscillators were fabricated 
as closely identical as possible, each attractor remains stable for a set of initial conditions and for a reasonably long 
time before drifting to another attractor. The individual attractor remains stable for a longer time and the drifting 
is slower when the mismatch is much smaller. On the contrary the drifting becomes faster for increasingly larger 
mismatch. This clearly indicates the existence of EM in a coupled system. Thus, given two closely identical systems, 
one can really observe various stable attractors for a change of initial conditions. On the other hand, this parameter 
mismatch could be taken care of from the beginning by giving a special attention to the design of coupling. By this 
way, one can obtain stable attractors avoiding the effect of drifting. For a demonstration, we consider two Rossler 
oscillators with non-identical parameters. The first oscillator has parameters a, b , c. 


±1 = -X2-X3, 

(5a) 

X 2 = X\ + aX2, 

(5b) 

±3 — b + x 3 (x 1 - c). 

(5c) 

The second oscillator has been detuned by A mismatch, 

m = -yi ~ 2 / 3 , 

(6a) 

V2 = yi + (a + A a )y 2 , 

(6b) 

2/3 = b + A b + 2 / 3 ( 2 /! -c-A c ). 

(6c) 


We again fix the parameters in the chaotic regime (a = 0.2, b = 0.2, c = 5.7) for both oscillators. 

To gain more insight into the parameter mismatch we analyze the general setup elaborated in section II. Let an 
n-dimensional system be governed by the equations x = F(x, /j,); F : R n —» R n ; where fi depicts all the parameters of 
the system. Another uncoupled system whose parameters are different from the first one is described by y = F(y, fi') 
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FIG. 3: (a) Maxima of x\ as a function of initial condition X 02 and mismatch A a , A a = —0.04; —0.02; 0; 0.02; 0.04. Other initial 
conditions are set to ®oi = 0.0; X 03 = —0.3; yoi = 3.1; yo 2 = —2.0; yo 3 = —0.1. (b) Maxima of x\ as a function of mismatch A a , 
when £02 is fixed at 2.0. 

and all the parameters of this system are denoted by yJ with fi' = y + A M . When coupling these two systems, we 
define the errors between the variables as e = x — ay where the description of a is the same as before. However, the 
error dynamics is governed by e = F(x, y) — ccF(y, y') + Ui — cm 2 . 

Our target is that the detuning of parameters cannot destroy the synchronization manifold as well as the complex 
behaviour of the coupled system, namely, the EM. Suppose we can construct an error dynamics in the following form 
e = (M - M') H ( x, y) . By designing the controllers in such a way that they depend on the mismatch, we are able to 
ensure EM for rather large intervals of mismatch. Moreover the deviation from the parameter of the first oscillator 
can be either positive or negative. 

The nonidentical Rossler oscillators have already been described by the equations [5ja)-(c) and (Hla)-(c). For sim¬ 
plicity we are dealing with an error dynamics for which we set all the an to 1, i.e., e; = Xi — yi,i = 1,2,3. One choice 
of the error dynamics is, e\ = £ 263 , e 2 = A a ei and £3 = —x 2 ei — ce 3 . One parameter (’a’) exists in the equation of x 2 , 
the other two (’b’ and ’c’) are connected with ± 3 . For simplicity we have introduced the mismatch only in parameter 
a. Therefore A& = A c = 0.0. 

We have taken a Lyapunov function V = \{e\ + e§) >0 for which V = eie'i + 6363 = — ce§ which ensures global 
stability for c > 0. Following the above mentioned condition, e 2 becomes constant C for any real values of A Q . The 
differences of the controllers can be obtained from the differences of [5^a)-[6Xa), E|b)-[ 6 |b), [5jc)-l6jc) following the error 
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dynamics, 


U11-U21 = -2/2 - 2/3 + x 2 + £3 + x 2 (x 3 - y 3 ), 
u 12 -u 22 = -ei - ae 2 + A a ei + A a 2/ 2 , 
ui 3 -u 23 = -x 3 (xi - c) + 223(2/1 - c) 

-c(x 3 - 2/3) - £2(3:1 - 2/i)■ 


(7) 

( 8 ) 


We set tin, U 22 and U 23 to zero. 

Then the coupled equations become 

±1 = —x 2 —x 3 , (9a) 

£2 = xi + ax 2 + ((A a - l)(xi - 2 / 1 ) - a(x 2 - y 2 ) + A a y 2 ), (9b) 

£3 = b+ (xi - c)x 3 + (-x 3 (xi - c) + 2 / 3(221 ~ c) 

~c(x 3 - 2/3 ) - £ 2 (£i - 2/i))> (9c) 

2/1 = -J /2 - 2/3 - (£ 2 ( 3 : 3 - 2 / 3)+£2 - 2/2+£ 3 - 2 / 3 ). (9d) 

2/2 = 2/1 + (« + A q ) 2 / 2 , (9e) 

2/3 = & + 2 / 3 ( 2 /!-c). (9f) 


We have plotted the maxima of x 2 in Fig. 0(a) showing the bifurcation like sequence obtained by varying the initial 
value X 02 from -4 to +4 for different values of the mismatch (A a = —0.04; —0.02; 0; 0.02, 0.04). Each sequence for a 
particular parameter mismatch A a in Fig. O a ) illustrates the fact that there exists a unique attractor (chaotic or 
limit cycle) for each value of X 02 • It is clearly seen that the EM is not destroyed by the parameter mismatch, positive 
or negative. However, the particular attractor found for a certain initial condition Xq 2 depends on the value of the 
mismatch as shown in Fig. [31(b). We mention that robustness with respect to parameter mismatch is established 
with this example in which a partial CS is achieved by the choice of an = 1. In general, the basic framework of the 
coupling configuration would, of course, also allow for other choice of a such as -2 to establish a partial AS state with 
an amplification of the attractor. This choice does not disturb the emergence of EM and its robustness as well. 


V. DISCUSSION 


We investigated the phenomenon of EM that reveals coexistence of infinitely many attractors in two coupled 
oscillators. EM is manifested in coupled n-dimensional systems with two basic properties: (1) partial synchronization 
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(CS or AS or MS) of n —1 components (state variables) while the n—th components keep a constant distance depending 
on the initial conditions, (2) emergence of a conserved quantity that characterizes the synchronization manifold in 
which the attractor lives. In the long-term limit the state space appears as foliated into infinitely many synchronization 
manifolds, each of them containing at least one attractor for each emergent constant determined by the set of initial 
conditions. As explained in detail in the introduction, we emphasize here again, that all attractors studied in this 
paper are in fact relative attractors, since they are attractors within the synchronization manifold in which they exist, 
but marginally stable in the direction perpendicular to the synchronization manifold. Our investigation is based on 
the derivation of controllers that allow an appropriate design of coupling leading to the state of EM. In this paper, 
this controller design, which has been introduced earlier in 20|, has been extended to achieve partial AS or MS in 


which some, say m, of the components are completely synchronized while n — m— 1 components are antisynchronized. 
This generalization of the coupling scheme which is based on Lyapunov function stability is also employed to amplify 
or attenuate the amplitude of the attractors of the coupled system. This property could be useful in an experimental 
set up leading to either very large or very small amplitude attractors. We have shown that the choice of controllers is 
wide and have given an analytical estimate of the possible number of different coupling schemes which could realize 
the desired extreme multistability. For each choice of the Lyapunov function and each possible choice of nullifying 
the differences between the controllers of the two coupled systems, we can get at least 2 n — 2 different designs of the 
controllers. We have shown that infinitely many options of coupling functions are available for realizing EM in two 
coupled systems. This offers a great flexibility, particularly, in view of an experimental design, since experimental 
setups might have some restrictions on a possible coupling design. Offering different coupling options always enhances 
the possibility of an experimental realization of the EM. 

Additionally, again with respect to experimental design, we addressed the robustness of EM to parameter mismatch 
since, in reality, oscillators will never be identical as usually assumed for theoretical investigations. It turns out that, 
with a proper design of the controllers, a larger mismatch is allowed without disturbing the basic conditions to sustain 
EM in coupled systems. 

The controller design has thus been shown to be very effective in realizing EM and very flexible. When EM emerges 
in two coupled systems, a conserved quantity appears in the long-term limit. This conserved quantity depends on 
the initial condition of the two coupled systems and determines the synchronization manifold in which a particular 
attractor lives. However, it is still an open question, how this controller design leads to simple or complex dependence 
of the conserved quantity on the initial conditions. To clarify this question could be subject of further research. 
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